Tuning of EAG K + channel inactivation: Molecular determinants 
of amplification by mutations and a small molecule 
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Ether-a-go-go (EAG) and EAG-related gene (ERG) K + channels are close homologues but differ markedly in their 
gating properties. ERG1 channels are characterized by rapid and extensive C-type inactivation, whereas mamma- 
lian EAG1 channels were previously considered noninactivating. Here, we show that human EAG1 channels 
exhibit an intrinsic voltage-dependent slow inactivation that is markedly enhanced in rate and extent by 1-10 pM 
3-nitro-Af(4-phenoxyphenyl) benzamide, or ICA1 05574 (ICA). This compound was previously reported to have 
the opposite effect on ERG1 channels, causing an increase in current magnitude by inhibition of C-type inactiva- 
tion. The voltage dependence of 2 pM ICA-induced inhibition of EAG1 current was half-maximal at — 73 mV, 
62 mV negative to the half-point for channel activation. This finding suggests that current inhibition by the drug 
is mediated by enhanced inactivation and not open-channel block, where the voltage half-points for current inhibi- 
tion and channel activation are predicted to overlap, as we demonstrate for clofilium and astemizole. The muta- 
tion Y464A in the S6 segment also induced inactivation of EAG1, with a time course and voltage dependence 
similar to that caused by 2 pM ICA. Several Markov models were investigated to describe gating effects induced by 
multiple concentrations of the drug and the Y464A mutation. Models with the smallest fit error required both 
closed- and open-state inactivation. Unlike typical C-type inactivation, the rate of Y464A- and ICA-induced inactiva- 
tion was not decreased by external tetraethylammonium or elevated [K + ] e . EAG1 channel inactivation introduced 
by Y464A was prevented by additional mutation of a nearby residue located in the S5 segment (F359A) or pore 
helix (L434A) , suggesting a tripartite molecular model where interactions between single residues in S5, S6, and 
the pore helix modulate inactivation of EAG1 channels. 



INTRODUCTION 

Ether-a-go-go (EAG) K + channels were first discovered 
in Drosophila melanogaster, where the behavioral mu- 
tant gag was associated with spontaneous repetitive fir- 
ing in motor axons and increased transmitter release 
(Ganetzky and Wu, 1983). Molecular cloning revealed 
that eag encoded a K + channel subunit (Warmke et al., 
1991), and voltage-clamp studies in Xenopus laevis oo- 
cytes determined that these channels conduct a de- 
layed rectifier K + current (Briiggemann et al., 1993). 
A human homologue (hEAG) of Drosophila eag a.nd a 
related gene (hERG, human EAG-related gene [ERG] ) 
were later discovered by screening a human hippo- 
campus cDNA library (Warmke and Ganetzky, 1994), 
and in turn, a Drosophila homologue of hERGl was 
identified (Titus et al., 1997). Further cloning efforts led 
to the identification of additional EAG (EAG1 and EAG2) 
and ERG (ERG1-ERG3) genes in mammals (Gutman 
etal., 2003). 
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The physiological roles of EAG and ERG channels 
have been extensively studied but remain incompletely 
characterized. In mammals, EAG1 channels are strongly 
expressed in several brain regions, including olfactory 
bulb, hippocampus, cerebellum, and brain stem (Martin 
et al., 2008), but the precise physiological roles of these 
channels are still not clear. Increased expression of 
EAG1 has been linked to certain cancers and tumor cell 
lines, and EAG1 inhibition by blockers, antibodies, and 
siRNA decreases the proliferation of tumor cell lines 
(Hemmerlein et al., 2006; Gomez-Varela et al., 2007; 
Pardo and Sttihmer, 2008). Drosophila ergl channels are 
encoded by the seizure locus, and mutations cause a 
muscle paralysis associated with hyperactivity in the 
flight motor pathway (Titus et al., 1997). In mammals, 
ERG1 channels are prominently expressed in the cen- 
tral nervous system with functions that may overlap with 
EAG channels. ERG1 channels are also expressed in the 
heart where their physiological and pathophysiologi- 
cal functions have been well characterized. Cardiac 
ERG1 channels conduct Zkt, the rapid delayed rectifier 



© 2012 Garg et al. This article is distributed under the terms of an Attribution-Noncom- 
mercial-Share Alike-No Mirror Sites license for the first six months after the publication 
date (see http://www.rupress.org/terms). After six months it is available under a Creative 
Commons License (Attribution-Noncommercial-Share Alike 3.0 Unported license, as de- 
scribed at http://creativecommons.Org/licenses/by-nc-sa/3.0/). 



The Rockefeller University Press $30.00 
J. Gen. Physiol. Vol. 140 No. 3 307-324 
www.jgp.org/cgi/doi/1 0. 1 085/jgp.201 21 0826 



307 



K + current (Sanguinetti et al., 1995; Trudeau et al., 

1995) that has a key role in determining duration of the 
plateau phase and the rate of repolarization of action 
potentials in cardiomyocytes. Loss of function muta- 
tions in human ERG1 prolong the action potential and 
cause long QT syndrome, a heart condition associated 
with potentially fatal cardiac arrhythmia (Curran et al., 
1995; Sanguinetti et al., 1996). 

Although closely related, the biophysical properties 
of EAG and ERG channels differ from one another. 
During a depolarizing voltage-clamp pulse, the rate of 
EAG current activation is relatively fast and exhibits a 
prominent Cole-Moore shift, i.e., delayed and slowed 
activation when currents are elicited from more nega- 
tive holding potentials (V^s) (Terlau et al., 1996; Lin 
and Papazian, 2007) . Drosophila eagl currents inactivate 
slowly (Brtiggemann et al., 1993; Robertson et al., 

1996) , whereas time-dependent inactivation was not de- 
tected for mouse EAG1 channels expressed in Xenopus 
oocytes (Robertson et al., 1996). ERG1 channels acti- 
vate more slowly than EAG1 and exhibit an exception- 
ally strong and rapid C-type inactivation (Smith et al., 
1996; Spector et al., 1996) . Inactivation gready reduces net 
repolarizing outward current at positive transmembrane 
potentials, permitting the prolonged plateau phase of 
cardiomyocyte action potentials typical for humans and 
other large mammals (Sanguinetti and Tristani-Firouzi, 
2006; Rosatietal., 2008). 

The mechanistic and structural basis of C-type in- 
activation has been explored in several K + channels 
using site-directed mutagenesis, x-ray crystallography, 
and molecular dynamics simulation. It has been pro- 
posed that C-type inactivation reduces K + conductance 
via a subtle voltage-dependent change in conformation 
of the selectivity filter that can be preceded by a transient 
loss of K + selectivity (Lopez-Barneo et al., 1993; Zhou 
et al., 2001; Cuello et al., 2010b). C-type inactivation is 
slowed by elevated [K + ] e and extracellular TEA (Grissmer 
and Cahalan, 1989; Choi et al., 1991; Lopez-Barneo 
et al., 1993), and mutations near the selectivity filter 
(e.g., S620T or S631A in human ERG type 1 [hERGl]) 
can greatly attenuate or eliminate this gating process 
(Schonherr and Heinemann, 1996; Suessbrich et al., 
1997; Zou et al., 1998). Recently, several synthetic com- 
pounds have been discovered that suppress C-type in- 
activation of hERGl channels and thereby cause a 
striking voltage-dependent increase in the magnitude 
of outward currents. One of these compounds, 3-nitro- 
jV-(4-phenoxyphenyl) benzamide, or ICA105574 (ICA), 
shifts the half-point of hERGl inactivation by +180 mV 
at a concentration of 2 uM (Gerlach et al., 2010), rival- 
ing or surpassing the effects of inactivation-perturbing 
point mutations. The effects of ICA on closely related 
EAG1 channels have not yet been determined. 

Here, we define the biophysical properties and struc- 
tural determinants of human EAG gene type 1 (hEAGl) 



channel inactivation using site-directed mutagenesis 
and pharmacological modulation of channels heterolo- 
gously expressed in Xenopus oocytes. We show that wild- 
type (WT) hEAGl channels exhibit an intrinsic slow 
inactivation that reduces peak outward current ampli- 
tude by <10%. Based on the effects of ICA on the re- 
lated hERGl channel described above, we assumed that 
ICA would cause a minor increase in hEAGl current 
magnitude consistent with removal of its subde inactiva- 
tion. Instead, we found that ICA inhibited hEAGl in a 
voltage-dependent manner consistent with an enhance- 
ment of its intrinsic slow inactivation gating process. 
While trying to understand the molecular mechanisms 
of these opposite pharmacological activities of ICA, we 
discovered that hEAGl inactivation was also markedly 
enhanced by nonconserved mutations of an aromatic 
residue (Tyr464) located in the S6 segment. Unlike 
typical C-type inactivation, the rate of Y464A- or ICA- 
induced inactivation was not slowed by elevated [K + ] e or 
external application of TEA. Moreover, in hERGl, the 
mutation of Tyr652 (equivalent to Tyr464 in hEAGl) 
does not appreciably alter channel inactivation. Markov 
models were developed that recapitulate the main fea- 
tures of intrinsic inactivation and its enhancement by 
the Y464A mutation or ICA. A homology model of the 
hEAGl channel pore domain was used to identify two 
residues in close proximity to Tyr464: Phe359 in S5 and 
Leu434 in the pore helix. Inactivation of Y464A EAG1 
channels was prevented by a second site mutation of 
either residue. Moreover, opposite to WT hEAGl chan- 
nels, but similar to WT hERGl, F359L hEAGl channels 
were activated by ICA. Collectively, these findings sug- 
gest that interactions between Tyr464, Phe359, and 
Leu434 are crucial molecular determinants of hEAGl 
channel inactivation. Similar to hEAGl channels, in- 
activation of KCNQ1 channels is also subtle (only 
clearly revealed by examination of tail current kinet- 
ics; Tristani-Firouzi and Sanguinetti, 1998), and point 
mutations of residues equivalent to Tyr464 and Leu434 
in hEAGl and a residue in S5 greatly accentuate 
KCNQ1 inactivation (Seebohm et al., 2005). Collec- 
tively, our present and previous site-directed mutagen- 
esis studies suggest a similar tripartite molecular model 
for slow inactivation of hEAGl and KCNQ1 channels, 
and that although hEAGl and hERGl channels are 
highly related, they do not share a common mecha- 
nism of inactivation. 

MATERIALS AND METHODS 

Molecular biology 

Human EAG1 (KCNH1; National Center for Biotechnology 
Information Reference Sequence: available from GenBank/ 
EMBL/DDBJ under accession no. NM_002238.3) cDNA cloned 
into psGEMHE oocyte expression vector was provided by the late 
Dr. Dennis Wray. Mutations in WT hEAGl cDNA were made 
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using the QuikChange site-directed mutagenesis kit (Agilent 
Technologies) and verified by DNA sequence analyses. Plasmids 
were linearized by NotI, and cRNA was prepared with the mMes- 
sage mMachine T7 kit (Life Technologies) . The concentration of 
cRNA was quantified by RiboGreen assay (Life Technologies) . 

Oocyte isolation and cRNA injection 

The procedure used to harvest oocytes from Xenopus was approved 
by the University of Utah Institutional Animal Care and Use 
Committee. To remove oocytes, frogs were anesthetized with 
0.2% tricaine methane sulfonate, and a small abdominal incision 
was made to remove ovarian lobes. The incision was sutured 
closed, and the frog was returned to its aquarium. A maximum of 
four surgical procedures with 1-mo intervals were performed on 
each frog. To isolate individual oocytes and remove their follicle 
cell layer, ovarian lobes were first manually dispersed with forceps, 
and then treated for 90-120 min with 1 mg/ml of type II collage- 
nase (Worthington Biochemical Corp.) in Ca 2+ -free ND96 solu- 
tion that contained (in mM): 96 NaCl, 2 KC1, 1 MgCl 2 , and 
5 HEPES, with pH adjusted to 7.6 with NaOH. Stage IV and V 
oocytes were injected with WT or mutant cRNA and incubated for 
1-5 d at 18°C in Barth's saline solution before use in voltage- 
clamp experiments. Barth's solution contained (in mM): 88 NaCl, 
2 KC1, 0.41 CaCl 2 , 0.33 Ca(N0 3 ) 2 , 1 MgS0 4 , 2.4 NaHC0 3 , 10 HEPES, 
1 pyruvate, and 50 mg/L gentamycin, with pH adjusted to 7.4 
with NaOH. 

Two-electrode voltage-clamp protocols and data analysis 

Ionic currents were recorded from single oocytes using a stan- 
dard two-microelectrode voltage-clamp technique (Goldin, 1991; 
Stuhmer, 1992). Oocytes were placed in a recording chamber 
(^0.4~ml volume) perfused at 2 ml min -1 with KCM211 solution 
at 22-24°C. KCM211 solution contained (in mM): 98 NaCl, 2 KC1, 
1 CaCl 2 , 1 MgCl,, and 5 HEPES, with pH adjusted to 7.6 with 
NaOH. A GeneClamp 500 amplifier, Digidata 1322A data acquisi- 
tion system, and pCLAMP 8.2 or 9.0 software (Molecular Devices) 
were used to produce command voltages and to record current 
and voltage signals. To reduce leakage of KC1 from micropipettes 
into the oocytes, agarose-cushion microelectrodes were fabricated 
as described previously (Schreibmayer et al., 1994). The resis- 
tance of microelectrodes after breaking their tip and sealing with 
1% agarose/3 M KC1 ranged from 0.2 to 0.6 Mil when backfilled 
with 3 M KC1. 

For experiments with TEA (Sigma-Aldrich) or when [KCl] e was 
increased, an equimolar concentration of NaCl was omitted to 
maintain constant osmolarity of the KCM211 solution. ICA (pur- 
chased from Sigma-Aldrich and AKos GmBH) was dissolved in 
dimethyl sulfoxide to prepare a 10-mM stock solution. Final [ICA] 
was obtained by dilution of the stock solution with KCM211 im- 
mediately before use for each experiment. 

During an initial equilibration period, oocytes were voltage 
clamped to a Vh of — 100 mV, and 1-s pulses to a test potential ( V t ) 
of 0 mV were applied every 10 s until current magnitude reached 
a steady-state level. During perfusion of the recording chamber 
with ICA solutions, the pulse interval was lengthened to 30 s. After 
currents achieved a new steady-state level in the presence of ICA, 
the relevant voltage-clamp protocols were repeated. 

To determine the voltage dependence of channel activation 
and inactivation, a dual-pulse protocol was used. Conditioning 
prepulses (V pre s) of 10-s duration to potentials ranging from 
— 130 to +20 mV were applied in 15-mV increments from a 14 of 
— 100 mV. After each prepulse, a 4.5-s test pulse was applied to 
+30 mV. Variations from this voltage-pulse protocol are de- 
scribed in Results and figure legends. Conductance (G) was cal- 
culated as the peak outward current in response to a prepulse 
(/p re ) divided by the corrected electrical driving force, as de- 
scribed in Eq. 1, 



p pre 

' ~ GHK[( V - E K ) ] ~ FV exp (F(V - E K ) / RT) =\ ' (1) 
RT exp(FV / RT)-1 

where E K is the reversal potential for K + flux (—96 mV in 2 mM 
of external K + ), Fis Faraday's constant, R is the gas constant, and 
Tis temperature in °K. The electrical driving force was corrected 
based on the Goldman-Hodgkin-Katz equation, defined as 
GHK[( V- E K )] (Clay, 2009), to account for unequal K + distribu- 
tion in the intracellular and extracellular solution. For each oo- 
cyte, G was plotted as a function of Vp re and fitted to a Boltzmann 
function. The fitted curve was extrapolated to more positive po- 
tentials to estimate maximum conductance ( G max ) . To determine 
the voltage dependence of current inactivation, the peak outward 
current during the test pulses to +30 mV (7 m ax, normalized to 1.0) 
was plotted as a function of Vp re , and the relationship was fitted 
with the Boltzmann function (Eq. 2): 

G _ 1 

G~ ~ 1 + exp (-zF(V pre - V 0jS )/RT)' 

where z is the effective valence, and is the potential at which the 
current is half-inactivated. In some experiments, adjusted / max was 
estimated by fitting the currents during each test pulse with an ex- 
ponential function and extrapolated to the onset of the test pulse. 

Offline data analysis was performed with Clampfit 9.0 (Molecular 
Devices), Origin 8.2 (OriginLab), and Excel (Microsoft) software. 
All data are expressed as mean ± SEM (n = number of oocytes) , and 
statistical significance was evaluated by two-way ANOVA and Tukey's 
HSD test, or the Student's t test where appropriate (P < 0.05 was 
considered a statistically significant difference). 

Biotinylation and Western blot 

The Pierce cell surface protein isolation kit (Thermo Fisher Sci- 
entific) was used for biotinylation of surface proteins in oocytes. 
Groups of oocytes (40-50) were injected with 20 ng of WT or 
mutant hEAGl cRNA and maintained for 4 d at 18°C in Barth's 
saline solution. Oocytes were then incubated with 1 mg/ml of 
solution (Sulfo-NHS-SS-Biotin in ND96) to label membrane pro- 
teins. After 1 h at 4°C, the reaction was stopped by adding quench- 
ing solution. Oocytes were washed twice with ND96 solution and 
once with Tris-buffered saline, and then incubated in 500 pi of 
lysis buffer with protease inhibitors for 15 min. After homogeniza- 
tion and incubation on ice for 15 min, the buffer volume was in- 
creased to 1 ml and lysate was centrifuged at 10,000 gfor 10 min. 
The clarified supernatant was dialyzed against NeutrAvidin Aga- 
rose beads (Thermo Fisher Scientific) for 3 h at 4°C. The beads 
were collected via centrifugation at 2,000 gfor 1 min and washed 
thoroughly with lysis buffer. After the final wash, samples were 
eluted by heating the beads to 60°C for 20 min in SDS-PAGE load- 
ing buffer. Samples were subjected to SDS-PAGE electrophoresis 
and analyzed by Western blotting with a 1:1,000 dilution of rabbit 
anti-hEAGl antibody (Abeam) and a 1:4,000 dilution of horserad- 
ish peroxidase-conjugated goat anti-rabbit secondary antibody 
(Santa Cruz Biotechnology, Inc.). Samples were also analyzed for 
calnexin (ER marker) with a 1:5,000 dilution of rabbit primary 
antibody (Abeam) and a rabbit secondary antibody (Santa Cruz 
Biotechnology, Inc.). GAPDH and G|3 were detected with a 
1:2,000 dilution of mouse primary antibody (Abeam) and a 1:500 
dilution of mouse antibody (Santa Cruz Biotechnology, Inc.), re- 
spectively, and a mouse secondary antibody (Jackson Immuno- 
Research Laboratories, Inc.). 

Molecular modeling 

The Swiss-Model Protein Structure Modeling Server at the Ex- 
PASy website was used in the automated mode (Arnold et al., 
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2006) to construct a homology model of the hEAGl channel based 
on the KcsA K + channel structure (Protein Data Bank accession no. 
1BL8) as a template. Structure diagrams were generated using 
PyMOL Molecular Graphics software (version 1.2.8; Schrodinger) . 

Markov modeling of EAG1 currents 

Mathematical models of EAG1 channel currents were developed 
to reconstruct experimental data. The explored models differ 
in their topology, as depicted in Figs. 10 and S9. Each model 
includes a Markovian description of channel states and the 
Goldman-Hodgkin-Katz current equation for K + current through 
a population of channels based on their states. In the Markov 
models, forward rates a and backward rates (3 were defined as 
dependent on the transmembrane voltage, V m : 



z a V m F/RT 



ipV m F/RT 



(3) 



(4) 



with the rates a 0 and p 0 at 0 mV, and the charges z„ and zp. Other 
transition rates (a C io, a CB , 0i O i, Pcio, Pcb, Poi, K y) had a 
constant value. Current conducted by EAG channels (/eag) was de- 
scribed by the Goldman-Hodgkin-Katz current equation for K + : 



/ =OP. 

EAG 



V n F[K + ],-[K + lexp(-V m F/RT) 
" G RT l-exp(-V m F/RT) 



(5) 



with permeability Peag, extracellular potassium concentration 
[K + ] e , intracellular potassium concentration [K + ]i, and the prob- 
ability that channels are in an open state O. 

Parameters of the Markov models were determined by numeri- 
cal fitting of feature vectors extracted from experimental and simu- 
lation data, f s and /„, respectively. The fitting was performed for 
experimental data from WT hEAGl at 0 (Control) , 2, 5, and 10 pM 
ICA, and from Y464A hEAGl. The fit error E was defined as: 

(6) 



II f ■ - f ■ II 1 

|| Jm,i Je,% llg 



fe,i 



+ (1 - Max £c/ + (Q Mra - Max £0/ , 



12 J 



with the number of features n, the Euclidean norm ||- • -|| 2 , the maxi- 
mal closed-state probability at the end of the prepulse (—90 mV) 

i 

the maximal open-state probability during the test pulse to +30 mV 
Max Z°/> 



current traces elicited by stepping to —130, —100... +20 mV. The 
resulting fitting parameters were applied for a final fitting using 
the complete dataset. In each of these fitting procedures, the 
number of iterations was eight. In each iteration, 480 perturbed 
parameter sets were evaluated. Perturbations were initially in the 
range of 50-150% of the original parameter value. After each it- 
eration, the perturbation range was decreased to a final range of 
98-102%. The 48 parameter sets with the smallest fit error were 
optimized using the steepest descent method (Press et al., 1992) 
and used for the parameter set perturbation in the next step. 

For the Control model, four degrees of freedom of the closed- 
inactivated and open-inactivated transitions were investigated. 
Rate coefficients oi CI0 and Pcio were kept identical to a 0 i and p OI , 
respectively. In addition, oi C n, a C i2, and a C B were identical as were 
Pen, Pcg. an <4 Pcb- For the ICA and Y464A models, linear interpo- 
lation was applied for a CI0 and p CI0 to a OI and p OI , respectively. 

Simulations with the models and their parameterization were 
performed using Matlab R201 lb (The Mathworks, Inc.) . The simu- 
lations were performed with a variable order method for solution 
of ordinary differential equations based on numerical differentia- 
tion formulas (Matlab function odel5s). An initial time step At of 
1 ms was chosen. Initial values of the state variables are listed 
in Table SI. For all models, the temperature was set at 293°K, 
[K + ] 0 = 2 mM and [K + ]j = 85 mM. All calculations were performed 
in double precision. Numerical fitting procedures were accelerated 
by using the Matlab Parallel Computing Toolbox. 

Online supplemental material 

Fig. SI illustrates the slow onset of open-state inactivation of WT 
hEAGl channels expressed in Xenopus oocytes. Fig. S2 shows that 
unlike typical Gtype inactivation, elevated [K + ] e does not signifi- 
cantly alter the rate of inactivation for Y464A hEAGl channel cur- 
rent. Fig. S3 compares the onset of fast inactivation for Y464A 
hEAGl channel current before and after treatment of an oocyte 
with 1 pM ICA. Fig. S4 shows that 10 mM of external TEA reduces 
current magnitude but does not slow the rate of inactivation for 
Y464A hEAGl channel current. Fig. S5 shows that 10 mM of exter- 
nal TEA does not slow the rate of ICA-induced inactivation of WT 
hEAGl channel current. Fig. S6 shows a Western blot for whole cell 
fraction of cRNA-injected and -uninjected oocytes. Fig. S7 shows 
that F359A and F359A/Y464A mutant hEAGl channels do not in- 
activate. Fig. S8 is a schematic representation of the five- and six- 
state Markov models for EAG1 channel gating. Fig. S9 presents 
simulated FAG1 currents using preliminary Markov models. Fig. 
S10 compares the rates of onset of activation and inactivation of 
WT hEAGl channel currents for oocytes treated with 2 and 10 pM 
ICA. Table SI lists features for model fitting and their fit errors. 
Table S2 presents the initial values for channel states used in differ- 
ent models. Table S3 lists the rate constant parameters for models 
of WT EAG1 channels (Control; 2, 5, and 10 pM ICA) and Y464A 
EAG1 channels. The online supplemental material is available at 
http://www.jgp.org/cgi/content/full/jgp.201210826/DCl. 



and a predefined maximal open-state probability 0 Max . Features 
were measured for /eag at different voltages and included maxima 
and parameters from exponential fits. The fit error E included 
functions of closed and open states to select parameter sets that 
cause specified open probabilities at activating voltages and high 
probability of closed states at low voltages, respectively. 

The fitting procedure was multi-scale and iterative, and it com- 
bined random perturbation and selection of parameter sets with 
the smallest fit error followed by their optimization. A multi-scale 
procedure was developed to establish estimates for initial values 
of the fitting parameters. A first estimate was determined by 
fitting to reduced data, i.e., features from current traces elicited 
by stepping to voltages of —130, —85, —40, and +5 mV. The 
estimates were used for a second fitting that used features from 



RESULTS 

Voltage-dependent inactivation of hEAGl channels 

To characterize the voltage dependence of hEAGl 
channel activation and inactivation, a double-pulse volt- 
age-clamp protocol was used. l^, re s of 10-s duration were 
applied in 1 5-mV increments to potentials ranging from 
— 130 to +20 mV. After each prepulse, a test pulse was 
applied to +30 mV to assess channel availability. Under 
control conditions, hEAGl currents activated with a V pTe ^ 
—40 mV did not exhibit any time-dependent decay in 
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magnitude (Fig. 1 A, top). However, the peak current 
measured at the test pulse of +30 mV (7 max ) decreased 
progressively as a function of T^ re , indicative of voltage- 
dependent inactivation. The extent of inactivation was 
small and only reduced 7 max elicited after a prepulse to 
+20 mV by 5-10%. The voltage dependence of hEAGl 
channel inactivation determined under control condi- 
tions is plotted in Fig. 1 B (closed squares) and fitted 
with a Boltzmann function to determine Vo.5 ( — 81.4 ± 
3.5 mV) and z (1.49 ± 0.32; n = 7). By comparison, the 
activation of hEAGl currents (Fig. 1 B, closed diamonds) 
developed over a more positive range of voltage, with a 
threshold near —55 mV, a of —4.6 ± 0.7 mV, and z of 
1.83 ± 0.06. Collectively, these findings indicate that 
hEAGl channels inactivate from both closed and open 
states. Control currents shown in Fig. 1 A do not exhibit 
a time-dependent decay expected for open-state inacti- 
vation, probably because the onset of inactivation is 
very slow and of minor proportion, and overlaps with a 
slow component of activation. A time-dependent decay 
of current, consistent with very slow open-state inacti- 
vation, was only clearly evident when the test pulse to 
+30 mVwas extended to 30 s (Fig. SI). 

Inhibition of hEAGl current by ICA 

ICA strongly inhibits inactivation and thereby in- 
creases the magnitude of hERGl channel currents, 
especially at positive potentials (Gerlach et al., 2010; 
Garg et al., 2011). We tested the effects of this com- 
pound on hEAGl, with the expectation that it might 
also inhibit intrinsic inactivation and increase cur- 
rent magnitude. Instead, we found that ICA appeared 
to enhance inactivation of hEAGl. As shown in the 
bottom of Fig. 1 A, after a 30-min incubation with 
2 uM ICA, peak currents were reduced in magnitude 
and exhibited further time-dependent decay during 
both the conditioning and test pulses. The prepulse 
voltage dependence of current inhibition by ICA was 
quantified by plotting 7 max as a function of Vp le . The 
relationship was fitted with a Boltzmann function 
(Fig. 1 B, open squares) to determine V0.5 ( — 66.0 ± 
4.1 mV) and z (1.15 ± 0.15; n=7). To account for the 
relatively rapid current inhibition compared with the 
slow onset of hEAGl activation at +30 mV, the decay 
phase of test-pulse currents was fitted to an exponen- 
tial function and an adjusted 7 max was estimated by 
extrapolation of the exponential fit to the onset of 
the test pulse. The adjusted 7 max was plotted as a func- 
tion of V pre (Fig. 1 B, open circles). With this adjust- 
ment, normalized 7 max after the conditioning pulse to 
-130 mV with 2 uM ICA matched predrug 7 max . The 
adjusted 7 ma x-Vpre relationship had a V0.5 of —73.3 ± 
3.9 mV and a z of 1.11 ±0.16. 

To determine the effects of ICA on hEAGl kinetics, 
currents recorded during the 10-s Vp re s were fitted with 
a mono- (control) or bi-exponential (ICA) function. 



ICA at 2 pM did not affect the rate of fast activation of 
hEAGl (Fig. 1 C, squares). The rate of ICA-induced cur- 
rent decay was voltage dependent, and the time constant 
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Figure 1 . Intrinsic inactivation of hEAGl channels is enhanced 
by ICA. (A) hEAGl currents recorded before (Control, middle) 
and after 2 pM ICA (bottom) using the voltage-pulse protocol 
shown in the top. V h was — 100 mV, and V prl . was 10 s in duration 
and ranged from — 130 to +20 mV, applied in 15-mV increments. 
After each V wc , a test pulse was applied to +30 mV. Zero current 
level is indicated by an arrow on the bottom left of the current 
traces in this panel and in all other figures. (B) Voltage depen- 
dence of hEAGl activation (control, ♦) and inactivation in the ab- 
sence (■) and presence of 2 pM ICA (□). For inactivation curves, 
normalized 7 max was plotted as a function of V prc . For the activation 
curve, /p r e normalized for electrical driving force was plotted as a 
function of Vp le . Data (ra = 7) were fitted with a Boltzmann func- 
tion (curves). For conuol inactivation, V 0 3 = —81.4 ± 3.5 mV and z = 
1.49 ± 0.32; for control activation, V os = -4.6 ± 0.7 mV and z = 1.83 ± 
0.06. For inactivation in the presence of ICA, V 0 3 = —66.0 ± 
4.1 mV and z = 1.15 ± 0.15. O, extrapolated 7 max in the presence 
of 2 pM ICA: Vo. 5 = -73.3 ± 3.9 mV and z= 1.11 ±0.16. (C) Plot of 
time constants for onset of activation (control, ■; 2 pM ICA, □) 
and inactivation (2 pM ICA, O) during the 10-s Vp re s (n = 8). 
(D) Voltage dependence of inactivation in the absence and pres- 
ence of three concentrations of ICA. The normalized / max -l^, re re- 
lationship for the test pulse is plotted for control (■), 2 pM ICA 
(solid curve; same as data plotted in B) , 5 pM ICA (O) , and 10 pM 
ICA (O). For reference, the voltage dependence of hEAGl activa- 
tion (B) is replotted as a dashed curve. For conuol inactivation, V0.5 = 
-82.6 ± 2.6 mV and z = 1.26 ± 0.36 (n = 16). For 5 pM ICA, 
V 0Ji = -80.0 ± 4.4 mV and z = 1.65 ± 0.27 (n = 5); for 10 pM ICA, 
V0.5 = -89.5 ± 4.7 mV and z = 1.39 ± 0.17 (n = 7). (E) Plot of the 
voltage dependence of hEAGl control activation and inactivation 
after 2 pM ICA. For inactivation, minimum and maximum 7 max 
was set to 0 and 1, respectively. 
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varied from 1.6 s at -40 mV to 0.9 s at +30 mV (Fig. 1 C, 
open circles). Higher concentrations of ICA (5 and 
10 pM) further reduced current magnitude and shifted 
the Vo.5 for current inhibition to more negative poten- 
tials (Fig. 1 D, open symbols). After the most negative 
prepulse to —130 mV, 7 max was reduced by 60% with 
5 uM and 83% with 10 pM ICA. 

Prepulse-dependent inhibition of current by 2 pM 
ICA developed over a voltage range that was mostly neg- 
ative to the voltage range required for channel activa- 
tion (Fig. 1 E) . The slopes of these relationships were 
similar; however, the V0.5 for channel activation was 
61 mV more positive than the V0.5 for current inhibi- 
tion. This difference suggests that ICA primarily in- 
hibits current by altering channel gating (enhanced 
inactivation) as opposed to a pore block mechanism 
where the voltage dependence of activation and inhibi- 
tion would be expected to fully overlap (Delpon et al., 
1997; Gonzalez et al., 2001; Lee et al., 2010). Clofilium 
(Gessner et al., 2004) and astemizole (Garcfa-Ferreiro 
et al., 2004) are two established open-channel blockers 
of hEAGl channels. Therefore, we next determined 
and compared the voltage dependence of activation 
and inhibition of hEAGl by these two drugs. 



Voltage dependence of open-channel block by clofilium 
and astemizole 

The effects of clofilium and astemizole on hEAGl cur- 
rents were examined at concentrations of 0.5 and 
2 pM, respectively. Currents were elicited with a two- 
pulse inactivation protocol similar to that used to 
characterize ICA, except that the maximum V pre was 
increased to +80 mV in the presence of drug. The cur- 
rent decay induced by astemizole during 10-s K re s 
(Fig. 2, A and B) superficially resembled the effects of 
ICA. However, as shown in Fig. 2 C, the V0.5 for chan- 
nel block by astemizole ( — 15.7 ± 3.3 mV; n = 5) was 
similar to the V0.5 of channel activation ( — 17.8 ± 
1.4 mV). The voltage dependence of block was inverted 
and its full range was normalized from 0 to 1 (Fig. 2 D) 
to more clearly illustrate that the voltage dependency 
of these two relationships overlaps one another. Clo- 
filium also inhibited hEAGl currents, but the onset of 
block was slower (Fig. 3, A and B) when compared 
with astemizole. In the presence of clofilium, V pie was 
increased to 30 s to account for the slower kinetics of 
block. The V0.5 for current inhibition by clofilium was 
nearly the same (shifted +6.5 mV) as the V0.5 of chan- 
nel activation (Fig. 3, C and D) . Thus, as predicted for 
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Figure 2. Voltage-dependent hEAGl channel block by astem- 
izole. (A) Control hEAGl currents (bottom) elicited with voltage- 
pulse protocol shown in the top. VJ, was — 100 mV, and V fIC was 10 s 
in duration and ranged from —130 to +35 mV, applied in 15-mV 
increments. (B) Currents recorded from the same oocyte after 
steady-state block by 2 uM astemizole. (C) Voltage dependence 
of hEAGl activation (♦) and inactivation (■) in control, and 
inhibition by 2 pM astemizole (□). For inactivation and inhibi- 
tion curves, normalized / max was plotted as a function of V* pre . For 
the activation curve, 7 pre normalized for electrical driving force 
was plotted as a function of V pre . Data (n = 5) were fitted with 
a Boltzmann function (curves). For control inactivation, V0.5 = 
-88.5 ± 2.4 mV and z = 2.84 ± 0.62; for control activation, V 0 3 = 
— 17.8 ± 1.4 mV and z = 1.84 ± 0.01. For inactivation in the pres- 
ence of astemizole, V0.5 = -15.7 ± 3.3 mV and z = 1.51 ± 0.07. 
(D) Comparison of the voltage dependence of hEAGl activation 
(dashed curve) and block by 2 pM astemizole (□). 
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Figure 3. Voltage dependence of hEAGl channel block by 
clofilium. (A) Control hEAGl currents (bottom) elicited with 
voltage-pulse protocol shown in the top. Vh was —100 mV, and 
Vp re was 10 s in duration and ranged from —130 to +35 mV, ap- 
plied in 15-mV increments. (B) Currents recorded from same 
oocyte after steady-state block by 0.5 pM clofilium. Note that the 
duration of Vp re was increased to 30 s. (C) Voltage dependence 
of hEAGl activation (♦) and inactivation (■) in control, and 
inhibition by 0.5 pM clofilium (□). For inactivation and inhibi- 
tion curves, normalized 7 max was plotted as a function of V pre . 
For the activation curve, 7 pre normalized for electrical driving 
force was plotted as a function of V pre . Data (n = 5) were fitted 
with a Boltzmann function (curves). For control inactivation, 
V0.5 = —86.7 ± 4.5 mV and z = 3.12 ± 1.21; for control activation, 
V 0 .5 = - 12.4 ± 0.8 mV and z = 2.00 ± 0.02. For inactivation in the 
presence of clofilium, V 0 5 = -5.2 ± 1.7 mV and z = 2.13 ± 0.13. 
(D) Comparison of the voltage dependence of hEAGl activation 
(dashed curve) and block by 0.5 pM clofilium (□). 
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Figure 4. Y464A hEAGl channels exhibit enhanced inactivation. 
(A) Voltage-pulse protocol and corresponding Y464A hEAGl cur- 
rents. (B) Test-pulse currents (during time indicated by a solid 
bar in A) are shown on an expanded time scale. Test currents elic- 
ited after V pre of-130, -85, -70, and -40 mV are labeled. A fast 
inactivating component of current is evident when V prc > — 70 mV. 

these open-channel blockers and unlike ICA, the 
voltage dependence of channel activation is a direct 
determinant of the relative block induced by astemi- 
zole or clofilium. 

Mutation of Tyr464 in S6 segment enhances hEAGl 
channel inactivation 

Mutation of Tyr652 in S6 of hERGl to Ala (Y652A) does 
not appreciably alter the kinetics of channel activation 
or inactivation (Fernandez et al., 2004) but greatly at- 
tenuates the ability of ICA to inhibit inactivation (Garg 
et al., 2011). Thus, we wondered if the corresponding 
mutation in hEAGl (Y464A) would also prevent or im- 
pact the ability of ICA to alter hEAGl channel gating. 
Unexpectedly, the Y464A mutation alone was found to 



induce prominent voltage-dependent inactivation of 
hEAGl (Fig. 4 A) . Moreover, inactivation was composed 
of multiple components that were revealed in a V pre -de- 
pendent manner. From a Vp re of —130 to —85 mV, cur- 
rent activated in response to a pulse to +30 mV exhibited 
fast and slow time components of inactivation, with 
time constants of 0.57 ± 0.02 s and 3.2 ± 0.3 s (Fig. S2; 
n= 5). From a V pre of —70 mV or more positive, an ad- 
ditional and even faster component of inactivation was 
evident (Fig. 4 B). When the Vp re was —85 mV or more 
negative, this fast component of inactivation was masked 
by the slower rate of activation. In addition to a gready 
enhanced inactivation, Y464A hEAGl channels ex- 
pressed poorly compared with WT channels, requiring 
injection of oocytes with ~14 ng cRNA to obtain equiva- 
lent current magnitudes after 1 d of incubation as that 
achieved with 0.2 ng cRNA for WT channels. 

In contrast to Y652A hERGl channels, which are in- 
sensitive to the activating effect of ICA (Garg et al., 
2011), Y464A hEAGl channels were even more sensi- 
tive to the inhibiting effect of ICA than WT channels, 
accelerating the rate of inactivation (Fig. S3) and caus- 
ing a greater tonic inhibition of current (Fig. 5, A and B) . 
For example, 1 uM ICA reduced 7 max ( V t = +30 mV and 
V h = -100 mV) by 35 ± 6% for WT channels (n = 5) 
versus 68 ± 4% (n = 8) for Y464A channels. 

The voltage dependence of activation and effects of 
ICA on inactivation {I mdL *-V pve relationship) of Y464A 
hEAGl are plotted in Fig. 5 C. The V0.5 for current acti- 
vation was —15.1 ± 2.1 mV (n = 8), a potential where 
inactivation was nearly complete. The I max -V pre relation- 
ship determined under control conditions (Fig. 5 C, 
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Figure 5. Effect of 1 uM ICA on Y464A 
hEAGl channel currents. (A) Control Y464A 
hEAGl currents elicited with two-pulse proto- 
col, where VJ, was — 100 mV and V pls was 10 s 
in duration and ranged from — 130 to +20 mV 
applied in 15-mV increments, each followed 
by a test pulse to +30 mV. (B) Currents re- 
corded after steady-state inhibition of hEAGl 
by 1 uM ICA. (C) Voltage dependence of ac- 
tivation (♦) and inactivation in the absence 
(■) and presence (O) of 1 uM ICA. For 
control activation, V 03 = —15.1 ± 2.1 mV and 
z = 1.31 ± 0.02 (n = 8). (D) Test-pulse cur- 
rents recorded at +30 mV after V plc to —1.30 
and —70 mV during control and after 1 uM 
ICA. Traces were fitted with two-exponential 
function and extrapolated back to the start of 
the pulse for control (blue curves) and after 
1 |iM ICA (red curves) . (E) Extrapolated / max - 
VJ, re relationships for control (■) and after 
1 uM ICA (□). Dashed curve is same activa- 
tion curve plotted in C. For control inactiva- 
tion, V0.5 = -73.3 ± 0.7 mV and z = 1.44 ± 0.05 
(n = 8). For 1 pM ICA, V 0£ = -94.3 ± 1.9 mV 
and z= 3.78 ±0.16 («=8). 



Garg et al. 313 



closed squares) and after 1 uM ICA (Fig. 5 C, open cir- 
cles) deviated from the normal sigmoidal shape. A 
probable cause of the reduced magnitude of I max at Vp re 
of less than —70 mV was the relatively slow onset of acti- 
vation of channels compared with their faster rate of 
inactivation. To account for this effect, the slower in- 
activating components of test-pulse currents were fitted 
to an exponential function, and an adjusted 7 max was es- 
timated by extrapolation of the exponential fit to the 
onset of the test pulse. This fitting procedure was per- 
formed for currents recorded before and after ICA 
(Fig. 5 D) . The voltage dependence of Y464A hEAGl 
inactivation (Fig. 5 E, closed squares) was determined 
by fitting the adjusted I mll x-V pTe relationship with a 
Boltzmann function ( V0.5 = -73.3 ± 0.7 mV and z = 1.44 
± 0.05). The relationship was nearly the same as that 
estimated for WT channels in the presence of 2 pM ICA 
(V0.5 = "73.3 ± 3.9 mV and z = 1.11 ± 0.16; Fig. 1 B). 
However, given the uncertainties inherent with exten- 
sive extrapolation required to estimate adjusted 7 max , 
especially in the presence of ICA (Fig. 5, D and E, open 
squares), the adjusted 7 m ax-Vpre relationship plotted in 
Fig. 5 E should be viewed with caution. Nevertheless, 
it is clear that ICA greatly enhanced inactivation of 
Y464A channels. 

Effects of external TEA and elevated [K + ] e 
on Y464A hEAGl inactivation 

The rate of C-type inactivation in Kv channels is slowed 
by high [K + ] e (Lopez-Barneo et al., 1993). In contrast, 
the rate of Y464A hEAGl channel inactivation was un- 
affected by elevation of [K + ] e from 2 to 20 mM, and 
slightly faster at 104 mM (Fig. S2). The rate of C-type 
inactivation is also slowed by external TEA (Choi et al., 
1991; Smith et al., 1996). In contrast, externally ap- 
plied 10 mM TEA did not slow inactivation of Y464A 
hEAGl, and in fact, at V t s between 0 and +40 mV, inac- 
tivation was accelerated by TEA (Fig. S4) . In addition, 
and similar to its effects on other Kv channels, 10 mM 
TEA caused a voltage-independent reduction of peak 
and steady-state current by ^50% (Fig. S4). TEA also 
had no apparent effect on the rate of ICA-induced 
inactivation of WT hEAGl channels (Fig. S5) . Col- 
lectively, these findings indicate that the prominent in- 
activation of hEAGl induced or enhanced by mutation 
(Y464A) or drug (ICA) is not mediated by a typical 
C-type mechanism. 

Molecular determinants of Y464A-induced inactivation 

To explore the structural basis of Y464A-induced inacti- 
vation in hEAGl, Tyr464 was substituted with other 
amino acids, including other aromatics (Phe and Trp) , 
a hydrophilic (Ser), or several hydrophobic (Leu, 
He, Met, and Val) residues. Preliminary experiments 
with these mutant channels indicated that some were 
nonfunctional (Y464I, Y464L, Y464M, and Y464W) or 



expressed very poorly (Y464V) . One possible explana- 
tion for the reduced or absent currents was that these 
mutant channels inactivated at much more negative 
potentials or a rate much faster than activation. It was 
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Figure 6. Mutation of Tyr464 in hEAGl has variable effects on 
inactivation. (A-E) WT or indicated mutant hEAGl currents re- 
corded from oocytes bathed in Mg 2+ -free external solution. WT 
currents were elicited with voltage-pulse protocol shown in the 
top: VJ, was —100 mV and V plc was 1 s in duration and ranged 
from —130 to +20 mV, applied in 15-mV increments. After each 
V prc , a 1.5-s test pulse was applied to +30 mV to measure channel 
availability. For some mutant channels V h and/or VJ,„ was varied 
as follows: Y464A, V h = —130 mV and V pre ranged from -160 to 
+20 mV; Y464S, V pre ranged from -130 to +95 mV; Y464V, V pK 
ranged from — 1 30 to +50 mV. (F) Voltage dependence of in- 
activation for mutant hEAGl channels (n = 4-6) as indicated. 
Boltzmann fits: Y464A, V 03 = -45.0 ± 1.3 mV and z = 1.14 ± 0.07; 
Y464S, V 03 =13.0 + 1.8 mV and z = 1.43 ±0.12; Y464V, V 03 = 
1 .6 ± 12.4 mV and z = 0.96 ± 0.13. (G; top) Western blots for hEAGl 
in the NeutrAvidin-captured cell surface protein fractions (mem- 
brane fraction) from a single batch of oocytes expressing WT or 
Tyr464 mutant channels. YL, Y464L; YI, Y464I; YM, Y464M; YV, 
Y464V) . The hEAGl monomer is represented by the 1 13-kD band; 
a higher band was also seen in all injected oocytes, but not in the 
uninjected oocytes (UN). Y464W (YW) channels and matched 
control (WT) channels were studied with a different batch of 
oocytes (less hEAGl protein; therefore, dimer band was absent) on 
a different gel. For all hEAGl constructs, each oocyte was injected 
with 20 ng cRNA. hEAGl antibody showed reactivity to multiple 
intracellular proteins in the whole cell fraction (WCF), overlap- 
ping the 11.3-kD hEAGl protein signal. Rightmost lane is for the 
WCF. (Bottom three panels) Western blots for calnexin, GAPDH, 
and G|3 from the same preparations as the top panel. Notice the 
absence of these three proteins in the membrane fraction. 
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previously reported that external Mg slows the rate of 
activation of EAG channels (Terlau et al., 1996) . There- 
fore, to accelerate channel activation and potentially 
unmask a fast component of inactivation, 1 mM Mg 2+ in 
the KCM2 1 1 extracellular solution was replaced by Ca 2+ . 
In addition, V pre was reduced to 1 s to allow us to pulse 
to more negative potentials without activating excessive 
endogenous currents. 

As illustrated in Fig. 6 A, WT channels activated more 
rapidly and displayed only minimal inactivation in 
the absence of external Mg 2+ . Similar to WT channels, 
Y464F hEAGl channel currents were robust and exhib- 
ited minimal inactivation (Fig. 6 B). For Y464A hEAGl , 
the faster rate of activation associated with the absence 
of Mg 2+ accentuated the fast component of inactivation 
(Fig. 6 C) . Y464S channels exhibited prominent inacti- 
vation, albeit requiring greater depolarization (Fig. 6 D) . 
Y464V hEAGl also required higher depolarization to 
elicit a modest inactivation (Fig. 6 E). The voltage de- 
pendence of inactivation of Y464F, Y464A, Y464S, and 
Y464V channels measured in the absence of external 
Mg 2+ is summarized in Fig. 6 F. 

All Tyr464 mutations other than Y464F exhibited 
poor or no functional expression. WT and Y464F 
hEAGl channel currents activated at +20 mV were typ- 
ically 5-7 uA when measured 24 h after injection of 
oocytes with 0.2 ng cRNA. In contrast, Y464A and 
Y464S mutants required larger amounts (14 ng) of 
cRNA to achieve similar levels of functional expres- 
sion. Nonaromatic and hydrophobic substitutions had 
the greatest effect on expression. Y464V channels ex- 
pressed at very low levels, and substitution with Trp, 
Leu, He, or Met resulted in nonfunctional channels re- 
gardless of incubation time (up to 6 d) when oocytes 
were injected with 14-20 ng cRNA. The Tyr464 muta- 
tions do not reduce or eliminate channel function be- 
cause of protein misfolding and/or reduced trafficking 
to the plasma membrane, as the surface expression of 
the Tyr464 mutant channels determined by a biotinyl- 
ation/Western blot assay was found to be similar to 
WT channels (Figs. 6 G and S6) . These findings sug- 
gest that Tyr464 mutants (other than Y464F) conduct 
poorly or not at all because they are highly or perma- 
nentiy inactivated. 

To further explore the molecular determinants of 
hEAGl inactivation, we mutated residues predicted to 
be in close proximity to Tyr464. A homology model of 
the hEAGl pore module was constructed based on the 
Streptomyces lividans K + channel KcsA crystal structure 
(Protein Data Bank accession no. 1BL8) . Tyr464 in 
hEAGl corresponds to Phel03 in KcsA, which was re- 
ported previously to interact with the pore helix to me- 
diate channel inactivation (Cuello et al., 2010a). The 
hEAGl model predicts that the side chain of Tyr464 
closely interacts with only two other residues: Leu434 in 
the pore helix and Phe359 in S5 of the same subunit 



(Fig. 7) . To test the importance of these two residues in 
channel inactivation, each was mutated singly to smaller 
or less hydrophobic residues, and in combination with 
Y464A to assay for intragenic rescue of gating. The con- 
servative mutation F359Y did not enhance inactivation 
nor did it rescue Y464A; i.e., F359Y did not prevent the 
Y464A-induced inactivation (Fig. 8 A). The noncon- 
served mutation F359L also did not enhance inactiva- 
tion but suppressed Y464A-induced inactivation; i.e., the 
double mutation F359L/Y464A restored normal gating 
(Fig. 8 B). F359Aalso rescued Y464A gating (Fig. S7). In 
addition, the K + selectivity of F359L/Y464A channels was 
reduced. After subtraction of linear leak, the reversal po- 
tential for these channels was —38.3 ± 1.4 mV (n = 7) in- 
stead of the usual —98 mV for WT or hEAGl channels 
with only a single mutation. Substitution of Leu434 with 
Met, Cys, or Ala had no effect on inactivation (Fig. 8, 
C-E, top). L434M/Y464A channels exhibited minimal 
but faster than normal inactivation, and L434C or L434A 
fully rescued the gating of Y464A channels (Fig. 8 C-E, 
bottom). Collectively, the site-directed mutagenesis find- 
ings described above suggest a tripartite model to explain 
the enhanced inactivation of Tyr464 mutant channels. 
We propose that in WT hEAGl , Tyr464 in S6 and Phe359 
in S5 strongly interact, presumably by ir-stacking. When 
this interaction is disrupted (e.g., by mutation of Tyr464 
to any residue other than Phe) , Phe359 can interact with 
Leu434 in the pore helix to strongly enhance intrinsic 
voltage- and time-dependent inactivation, presumably by 
a conformational change in the selectivity filter that re- 
duces ion conductance that is similar but not identical to 
Gtype inactivation. Substitution of Tyr464 with other 
hydrophobic residues (Leu, He, Met, and Val) has an 
even greater disruptive effect on the selectivity filter, per- 
haps via direct interaction with Leu434 in the pore helix 
in addition to release of Phe359 from the braking action 
exerted through the putative ir-stacking interaction 
with Tyr464. The Tr-stacking hypothesis predicts that 




Figure 7. Homology model of the pore domain for a single 
hEAGl subunit. Model used the crystal structure of KcsA (Pro- 
tein Data Bank accession no. 1BL8) as a template. Lateral view 
(A) and close-up view (B) of a subunit with Leu434 (base of pore 
helix), Tyr464 (S6), and Phe359 (S5) residues shown in space fill. 
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Figure 8. Intragenic rescue of Y464A-enhanced 
hEAGl inactivation by a second site mutation in 
either Phe359 (in S5) or Leu434 (in pore helix). 
(A-E) Currents for indicated mutant hEAGl chan- 
nels recorded during two-pulse inactivation pro- 
tocol: K = —100 mV, and V prc (10 s) ranged from 
— 130 to +20 mV, applied in 15-mV increments. 
A test pulse to +30 mV was applied after each V plc . 
The top set of traces in each panel is for channels 
with indicated single mutation; the bottom set of 
traces shows currents for same mutation introduced 
into the Y464A channel. All mutations except the 
conserved F359Y (A, bottom) and L434M (C, bot- 
tom) rescued inactivation of Y464A. Inward currents 
in some traces represent unsubtracted leak currents. 
Based on reversal potential of time-dependent cur- 
rents, only F359L/Y464A channels had reduced 
K + selectivity. Currents shown are representative of 
four to six oocytes for each mutant channel. 



substitution of Tyr464 with Tip, another aromatic resi- 
due, would not alter channel gating. Instead, we found 
that Y464W mutant channels trafficked normally to the 
cell surface but were nonfunctional. Tip is more bulky 
than Tyr or Phe, and this may prevent its TT-stacking with 
Phe359 and/or enhance direct interaction with Leu434 
to inactivate the channel similar to that proposed for the 
other nonaromatic hydrophobic substitutions. 

F359L hEAGl is activated by ICA 

As noted above, the voltage dependence of Y464A 
hEAGl inactivation is similar to that induced in WT 
channels by 2 uM ICA, suggesting a common under- 
lying mechanism. To explore this possibility, we deter- 
mined the effect of ICA on F359L mutant channels. 
F359L is one of the mutations that rescued gating of 
Y464A channels (Fig. 8 B) , and the equivalent mutation 
in hERGl (F557L) eliminated channel response to 30 uM 
ICA, presumably because it disrupted the drug-binding 
site (Garg et al., 2011). Based on its effects in hERGl, 
we assumed that F359L hEAGl channels would also be 
insensitive to ICA. Instead, we found that ICA increased 
F359L hEAGl channel currents by 5.4 ± 0.6-fold at 
10 pM (n = 3) without any evident change in kinetics 
(Fig. 9 A) . Thus, in contrast to F557L hERGl and op- 
posite to the effect on WT hEAGl, ICA is an agonist of 
F359L hEAGl channels. The mechanism of this activity 



will require further study but presumably results from 
an increase in channel open probability, perhaps by 
preventing F359L hEAGl channel transition into an 
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Figure 9. ICA activates F359L hEAGl channels and inhibits 
F359L/Y464A hEAGl channels. (A) F359L hEAGl currents re- 
corded under control conditions and after the application of 
10 uM ICA. (B) F359L/Y464A hEAGl currents recorded under 
control conditions and after the application of 10 pM ICA. In all 
panels, currents were elicited with two-pulse inactivation proto- 
col: V h = -100 mV, and (10 s) ranged from -130 to +20 mV, 
applied in 15-mV increments. A test pulse to +30 mV was applied 
after each V me . 



316 



EAG1 channel inactivation 




«CI 0 || PciO a Cli||Pcil Otc^jj Pci2 «C13 ||pCt3 a OI || Pol «Oi||Pc 



B 




a cio I I Pao a< || Pr : "c^jj Pci2 « ||( a o, | | P 




Figure 10. Schematics of (A) 12- and 
(B) 10-state Markov models of EAG1 chan- 
nel current. 



absorbing inactivated state. Opposite to channels har- 
boring only the F359L mutation, F359L/Y464A hEAGl 
channel currents were inhibited by 80.4 ± 0.02% (n= 7) 
with 10 uM ICA (Fig. 9 B). However, unlike WT chan- 
nels, F359L/Y464A hEAGl channel currents recorded 
in the presence of ICA did not exhibit time-dependent 
decay. This finding suggests that either ICA blocks the 
pore of these channels or primarily enhances closed- 
state inactivation. Thus, ICA is an activator of F359L 
hEAGl and an inhibitor of F359L/Y464A hEAGl chan- 
nels, indicating that Phe359 in S5 is required for ICA to 
induce open-state inactivation, as well as demonstrating 
that Tyr464 in S6 is required for the activator effect 
of ICA on F359L hEAGl channels. The latter finding is 
of particular interest because we reported previously 



that the equivalent Tyr in hERGl (Tyr652) is also 
required for ICA-mediated channel activation (Garg 
et al., 2011). Collectively, our present findings with 
hEAGl and previous findings with hERGl indicate that 
a strategically positioned Phe in S5 and a Tyr in S6 are 
critical determinants of ICA-induced alterations of 
channel inactivation. 

Markov modeling of hEAGl channel gating 

Mathematical models of WT EAG1 channel currents 
were developed to reconstruct experimental data. The 
first and simplest model considered four closed states 
and a single open state (Fig. S8 A) . This model was able 
to satisfactorily reproduce I-V relationships and current 
kinetics during the prepulses (Fig. S9 A); however, it 
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Figure 11. Comparison of measured and modeled currents for WT hEAGl. (A-D) Currents recorded from five to eight oocytes for 
each condition (Control [A] or ICA at 2 [B] ,5 [C] , and 1 0 uM [D] ) were averaged and normalized to their corresponding peak control 
current. Voltage-pulse protocol: V h = —90 mV and V^ re (10 s) ranged from —130 to +20 mV, applied in 30-mV increments; at 10 s, a test 
pulse to +30 mV was applied. (E-H) Simulated EAG currents (/) for control (E, 12-state model) and in the presence of 2 (F), 5 (G), and 
10 uM (H) ICA (all 10-state models). 
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did not reproduce the experimentally observed In- 
dependent differences in 7 max for test pulses to +30 mV. 
A six-state model that included a single inactivated state 
coupled to the single open state (Fig. S8 A) was also un- 
able to reconstruct differences in 7 max (Fig. S9 B). We 
next considered 10- and 12-state models that paired an 
inactivated state to each closed and open state. The 
12-state model featured 2 open states (Fig. 10 A); the 
10-state model had a single open state (Fig. 10 B) . Only 
the 12-state model (Fig. 11 E) accurately reproduced 
the slower phase of activation observed for WT hEAGl 
under control conditions (Fig. 1 1 A) . This difference is 
more obvious by inspection of test currents (at +30 mV) 
in Fig. S9 C (10-state model) and Fig. 11 E (12-state 
model) and reflected in the fit errors (Table 1). For the 
12-state model, the errors for individual features and 
the rate coefficients are presented in Tables S2 and S3, 
respectively. The modeled magnitudes of peak and end 
of pulse currents under control conditions are plotted 
in Fig. 12 (A and B) for the prepulses and in Fig. 13 (A 
and B) for the test pulses and closely match the experi- 
mentally measured values. An important feature of the 
Control model was accentuated forward rate coeffi- 
cients QL Cn , ot CI2 , and ctci3 versus forward rate coefficients 
a C io and a ol . 

Multiple models were considered for simulating cur- 
rents measured in the presence of 2, 5, and 10 pM ICA. 
Parameterization of the five-state model, lacking an in- 
activated state, and the six-state model, with an inacti- 
vated state, was not successful (Fig. S9, D and E). We 
next attempted to vary only C to I and O to I transitions 
but keeping ot 0 , z a, Po, and zp fixed as in the 10- and 
12-state control models; however, the resulting fit error 



TABLE 1 

Fit error E for the investigated models 

Model E E E E 

5-State 6-State 10-State 12-State 

Control 3.894 3.892 0.421 0.271 s 

2uMICA 230.212 0.696 0.430" 0.463 

5uMICA 981.122 1.965 1.472" 1.512 

lOuMICA 47.436 1.454 1.133 a 1.154 

Y464A 39.876 0.433 0.232 a 0.254 



a The smallest E for control and different ICA concentrations. 

was large. The fitting was improved by using the com- 
prehensive parameterization approach as for the con- 
trol model, and a linear interpolation of rate coefficients 
as described in Materials and methods. Results of pa- 
rameterization of the 10- and 12-state models are pre- 
sented in Figs. 1 1 D and S8 F, respectively. The fit error 
was similar for the 10- and 12-state models, and both 
of these were smaller than the 6-state model (Tables 1 
and S2). Therefore, the 10-state model was selected 
for further analysis of currents recorded in the pres- 
ence of ICA. 

The experimental data for 2, 5, and 10 pM ICA are 
presented in Fig. 11 (B-D) and are compared with 
model currents in Fig. 11 (F-H). The models repro- 
duce the [ICA] -dependent decrease in relative current 
magnitude, a manifestation of enhanced inactivation. 
For the three ICA models, the magnitudes of peak and 
end of pulse currents are plotted in Fig. 12 (C-H for 
prepulses) and in Fig. 13 (C-H for test pulses) and ad- 
equately match the experimentally measured values. 
In the three ICA models, the major effects of drug 




V(mV) V(mV) V (mV) V (mV) 

Figure 1 2. Features of measured and modeled WT hEAGl current during prepulses. (A-H) Plots of normalized maximum (Ipremax) and 
end (Ipr ee „d) currents during the prepulse for control (A and B) in the presence of 2 (C and D) , 5 (E and F) , and 10 pM ICA (G and H) . 
In each panel, experimentally measured currents (curves) are compared with modeled currents (X), all normalized to their respective 
-Tmax (peak current at V t of +30 mV) . 
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Figure 1 3. Features of measured and modeled WT hEAGl current during test pulses to +30 mV. (A-H) Plots of normalized maximum 
(Imax) and end (I m d) currents during test pulse for control (A and B) in the presence of 2 (C and D), 5 (E and F), and 10 pM ICA (G 
and H). In each panel, experimentally measured currents (curves) are compared with modeled currents (X), all normalized to their 
respective 7 max . 



compared with the Control model were faster forward 
rates for closed-state inactivation (a C n, oici2> an d a CIS ) 
and the forward rate a 0 i for open-state inactivation. For 
example, cice was increased from 0.0256 in Control to 
0.257 for 2 pM ICA, whereas the backward rate (Pce) 
was only increased from 0.016 to 0.089 (Table S3) . 

The 10- and 12-state models were equally effective 
(similar fit errors; Table 1) in reproducing Y464A 
hEAGl channel currents. Individual fit errors for the 
Y464A model are presented in Table SI. The experi- 
mentally measured Y464A currents and modeled cur- 
rents are compared in Fig. 14 (A and B). The model 
reproduces the mutation-induced change in peak cur- 
rents and enhanced inactivation (Fig. 14, C-F) but 
does not account for the very fast phase of inactivation 
measured at +30 mV after V pre > -70 mV (Figs. 4 B 
and 6 C). The major difference between the Y464A 
model and all other models was the slowed activation 
and deactivation rates (Table S3). In addition, com- 
pared with the Control model, the Y464A model had 
faster forward and backward rates for closed- and 
open-state inactivation. For example, ctci2 was increased 
from 0.0256 (Control) to 0.168 (Y464A), and the back- 
ward rate (Pgr) was increased from 0.016 (Control) to 
0.247 (Y464A) (Table S3). 

DISCUSSION 

Inactivation of EAG1 channels appears to be species 
dependent. Drosophila eagl channels inactivate sub- 
stantially (Briiggemann et al., 1993; Robertson et al., 
1996) , whereas inactivation is absent in mouse EAG1 
(Robertson et al., 1996) . Here, we report that human 
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Figure 14. Comparison of measured and modeled currents 
of Y464A hEAGl. (A) Currents recorded from 10 oocytes 
were averaged and normalized to I maK . Voltage-pulse protocol: 
Vt=-90mVand V pre (10 s) ranged from -130 to +20 mV, applied 
in 30-mV increments; at 10 s, a test pulse to +30 mV was applied. 
(B) Simulated EAG1 currents (I) for Y464A hEAGl. Features of 
the currents were extracted during the prepulse (C and D) and 
test pulse (E and F) . In each panel, experimentally measured cur- 
rents (curves) are compared with model currents (X), all normal- 
ized to their respective / max . 
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EAG1 channels exhibit a modest (5-10% reduction 
in current at +30 mV) and very slow voltage-depen- 
dent inactivation. Experimental data demonstrate 
and Markov modeling corroborates that hEAGl chan- 
nels inactivate from both closed and open states. 
Intrinsic inactivation was markedly enhanced by ICA, 
a previously described activator of hERGl channels. 
Intrinsic inactivation was also enhanced by mutation 
of Tyr464, equivalent to a key aromatic residue previ- 
ously implicated in KcsA channel inactivation (Cuello 
et al., 2010b). The molecular determinants for en- 
hanced inactivation induced by ICA or the Y464A mu- 
tation may be similar and are tentatively proposed to 
involve interaction between three residues located 
in close proximity to one another in S5 (Phe359) , S6 
(Tyr464) , and the pore helix (Leu434) in a single 
hEAGl subunit. 

Comparison of hEAGl inactivation with other K + channels 
Several modes of inactivation have been defined in 
K + channels, including N-, C-, P-, and U-type inactiva- 
tion (Kurata and Fedida, 2006). N-type inactivation 
was first identified in Shaker channels and is mediated 
by physical occlusion of the pore by the cytoplasmic 
N-terminal peptide chain (Hoshi et al., 1990). After 
truncation of the N terminus in Shaker, a much slower 
C-type inactivation was revealed (Hoshi et al., 1991). 
C-type inactivation is mediated by a conformational 
change in the selectivity filter that prevents or reduces 
K + conductance. A similar mechanism likely mediates 
the more rapid P-type inactivation (De Biasi et al., 
1993), and thus, the lumped term "P/ C-type" inactiva- 
tion is often used (Chen et al., 2000). A hallmark 
of C-type inactivation is that its rate is slowed by extra- 
cellular TEA or an increase in [K + ] e (Grissmer and 
Cahalan, 1989; Choi et al., 1991; Lopez-Barneo et al., 
1993; Smith et al., 1996). In contrast, external TEA in- 
creases currents in channels that exhibit P-type inacti- 
vation (De Biasi et al., 1993). ERG1 channels show 
prominent C-type inactivation, which is sensitive to ex- 
ternal TEA (Smith et al., 1996), [K + ] e (Wang et al., 
1997) , and to mutations of residues near the selectivity 
filter (Schonherr and Heinemann, 1996; Smith et al., 
1996; Suessbrich et al., 1997; Zou et al., 1998). In con- 
trast, we found that inactivation of the related hEAGl 
was not slowed by external TEA or elevated [K + ] e . 
Other types of inactivation (e.g., U-type) have been de- 
scribed in Kv2.1, Kv3.1, Kv4, and Shaker (Jerng and 
Covarrubias, 1997; Klemic et al., 1998, 2001; Kurata 
and Fedida, 2006), where channels inactivate from 
the preopen and open states and are insensitive, or in 
some cases their rate of onset is enhanced, by extracel- 
lular TEA or elevated [K + ] e . Further study is required 
to understand how slow inactivation in hEAGl differs 
or is similar to specific properties of inactivation de- 
scribed for other Kv channels. 



Markov models of EAG 

The Markov models presented here are the first quanti- 
tative description of EAG channel gating and currents. 
Four different model topologies were investigated using 
an automated stochastic optimization approach de- 
scribed previously for our study of hERGl gating cur- 
rents (Abbruzzese et al., 2010) . This approach evaluates 
a far greater number of parameterizations than are nor- 
mally investigated using the more common methodol- 
ogy of manual parameter adjustment. The approach 
allowed us to demonstrate that simple 5- and 6-state 
models were inadequate and that a 10- or 12-state model 
was more appropriate for description of EAG1 chan- 
nels. The 10- and 12-state models featured an inacti- 
vation state paired with each closed and open state 
because there was no compelling reason to exclude the 
likelihood that channels could inactivate from any 
particular state. Lacking experimental evidence to the 
contrary and to simplify parameterization, the rate co- 
efficients between inactivated states were identical to 
those of their paired closed or open states. For the con- 
trol model, the best fit was obtained using two open 
states to account for the slow phase of activation of 
currents elicited from the more negative V pre . For 
ICA or Y464A models, a single open-state model was 
as effective as a two open-state model because of the 
dominance of inactivation gating. There are several 
limitations of the models. First, some important experi- 
mental data are lacking, including single-channel con- 
ductance and kinetics, gating currents, and kinetics of 
drug binding. The lack of these measures limited our 
ability to further constrain the model and motivated 
our simplistic definition of parameters. Second, the 
analysis was constrained to four topologies; many others 
are possible but were not considered. 

ICA enhances inactivation of hEAGl channels 
The time- and voltage-dependent inhibition of hEAGl 
channel currents by ICA could result from either a di- 
rect consequence of binding to the central cavity (pore 
block) , or indirectly by enhancing the rate and extent 
of intrinsic inactivation. Astemizole (Garcfa-Ferreiro 
etal., 2004; Gessner etal., 2004) and clofilium (Gessner 
et al., 2004) have been described as direct pore blockers 
of EAG1 channels, whereas the phenylquinazolinone 
HMJ-53A (Chao et al., 2008) and the natural alkaloid 
rhynchophylline (Chou etal., 2009) inhibit delayed rec- 
tifier K + channels in neuroblastoma N2A cells by en- 
hancing C-type inactivation. We propose that ICA is a 
gating modifier and not a pore blocker for several rea- 
sons. First, the voltage dependence of hEAGl channel 
activation does not overlap with the voltage depen- 
dence of ICA-induced inhibition, as would be expected 
for an open state-dependent pore blocker. For ICA, the 
Vo.5 for inhibition of current was 62 mV more negative 
than the Vo.5 for current activation. Moreover, as relative 
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channel activation was increased from 20 to 90%, inhi- 
bition was only increased from 77 to 79%, inconsistent 
with open-channel block. In contrast, for two known 
pore blockers, astemizole (Garcfa-Ferreiro et al., 2004) 
and clofilium (Gessner et al., 2004), we found that the 
voltage dependence for channel activation and pre- 
pulse-dependent inhibition were the same. Second, the 
rate of current decay (inactivation rate) was only mod- 
estly accelerated when [ICA] was increased from 2 to 
10 uM, despite a great increase in tonic inhibition 
(Fig. S10). This finding is in contrast to the dramatic 
[astemizolej-dependent increased rate of time-depen- 
dent block of hEAGl (Garcfa-Ferreiro et al., 2004) chan- 
nels. Enhanced tonic inhibition caused by higher [ICA] 
is attributed to the trapping of channels in an absorbing 
inactivation state. Third, IG\ exerts its activator activity in 
hERGl by binding to a site outside of the central cavity, 
specifically to a hydrophobic-binding pocket located be- 
tween S5 of one subunit and S6 of a neighboring subunit 
(Garg et al., 201 1) . Because of the high sequence homol- 
ogy between the pore domains of hEAGl and hERGl, it 
is likely that ICA binds to a similar site in both channels. 
A key residue of the activator-binding site in hERGl is 
Phe557; F557L channels were insensitive to ICA (Garg 
et al., 2011). The homologous mutation in hEAGl (F359L) 
also dramatically altered drug action; instead of inhibi- 
tion, ICA activated F359L hEAGl. This key Phe residue is 
located in S5 away from the central K + conduction path- 
way. The side chain of Tyr464 in S6 faces the central cav- 
ity. Mutation of the corresponding residue in hERGl 
(Y652A) suppresses activity of many known pore block- 
ers, presumably because it is a key component of the 
binding site for these drugs (Mitcheson et al., 2000; 
Fernandez et al., 2004). In contrast, the mutation Y464A 
in hEAGl enhanced rather than reduced current inhibi- 
tion by ICA. Collectively, our present and previous find- 
ings strongly suggest that ICA does not bind to the central 
cavity of hEAGl. Instead, we propose that ICA modifies 
inactivation gating of hEAGl by binding to the same or 
similar hydrophobic pocket associated with ICA-medi- 
ated suppression of inactivation in hERGl channels. 

Mechanism of intrinsic hEAGl inactivation and its 
enhancement by Y464A or ICA 

Slow inactivation of WT hEAGl channels is voltage 
dependent and results in <10% reduction of current 
magnitude at positive potentials. It is difficult to ascertain 
whether the strong inactivation of Y464A hEAGl chan- 
nels is simply an exaggerated form of the intrinsic gating 
described for WT channels or is an additional and dis- 
tinct mutation-induced form of inactivation. However, it 
is noted that the voltage dependence and rate of inactiva- 
tion of Y464A hEAGl appear superficially similar to the 
intrinsic gating of Drosophila eagl. Perhaps coincidendy, 
the voltage dependence of inactivation of Y464A hEAGl 
is also similar to that induced by 2 uM ICA in WT hEAGl . 



One difference was that Y464A hEAGl inactivation in- 
cluded a very rapid component that was most evident 
when currents were elicited from a 14 of —70 mV or 
more positive. This fast component was not detectable 
with ICA-induced inactivation. Nonetheless, similar to 
ICA-induced inactivation, Y464A hEAGl inactivation 
was unaffected (or slighdy accelerated) by external TEA 
or high [K + ] e , implying a mechanism (s) different from 
C-type inactivation. 

Y464A hEAGl was more sensitive to the inhibitory ef- 
fects of ICA but suppressed the ability of ICA to activate 
F359L hEAGl channels. Likewise, the equivalent muta- 
tion in hERGl (Y652A) prevented channel activation by 
ICA. Together, these findings suggest that this key Tyr 
residue may be an allosteric mediator of both the activa- 
tor and the inhibitory effects of ICA on these related 
channels. The functional importance of Tyr464 can 
also be inferred from the functional consequences of 
mutation to residues other than Ala or Phe. Substitu- 
tion of Tyr with Ser or Val enhanced inactivation but 
also shifted the voltage dependence to more positive 
potentials. Substitution with the hydrophobic residues 
Leu, He, or Met generated nonfunctional channels that 
were nonetheless trafficked normally to cell surface, 
suggesting the possibility that Y464L/I/M channels are 
constitutively inactivated. Mutation of the equivalent 
S6 residue (Phe344) in KCNQ3 channels also dramati- 
cally decreased currents while having little or no effect 
on surface membrane expression of channel proteins 
(Choveau et al., 2012a). Homology modeling of these 
and other mutant KCNQ3 channels was used to con- 
clude that specific interactions between the pore helix 
and S6 can alter the structure of the selectivity filter to 
either favor or prevent K + conductance (Choveau et al., 
2012a, b). Thus, as discussed further below, an inter- 
action between a specific residue in S6 and the pore 
helix appears to mediate inactivation in KCNQ and 
other K + channels not directly related to EAGL 

Mutation of the residue equivalent to Tyr464 can af- 
fect inactivation gating in other Kv channels. Based on 
sequence alignments, Tyr464 in hEAGl is homologous 
to Phe 103 in KcsA, Ile470 in Shaker, and Tyr652 in 
hERGl channels. Phel03 in KcsA is proposed to be 
involved in bidirectional coupling between the activa- 
tion gate (S6 bundle crossing) and the inactivation gate 
(selectivity filter) (Cuello etal., 2010a). Mutation of Ile470 
in Shaker or Phe 103 in KcsA to smaller residues sup- 
presses inactivation (Holmgren etal., 1997; Melishchuk 
and Armstrong, 2001; Cuello et al., 2010a). In Shaker, 
substitution of Ile470 with a Tyr induces a complex in- 
activation with several components that show variable 
sensitivity to [K + ] e and external TEA (Klement et al., 
2008), and I470A channels have reduced inactivation 
(Cuello et al., 2010a). In KcsA, elegant structural stud- 
ies (Cuello et al., 2010a) have shown that movement 
of the inner (TM2) helix associated with activation is 
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linked with a tilting of Phel03 toward the pore helix, 
and via a network of H bonds among Trp67, Glu7l , and 
Asp80 causes inactivation by destabilization of the selec- 
tivity filter. However, opposite to the effects of the muta- 
tions Y464A in hEAGl or F344A in KCNQ3 (Choveau 
et al., 2012a), F103A KcsA has a reduced rate and ex- 
tent of inactivation compared with WT channels ( Cuello 
et al., 2010a). Finally, in contrast to Shaker, KcsA, 
KCNQ3, and hEAGl, the equivalent S6 mutation in 
hERGl (Y652A) has minimal effect on inactivation gat- 
ing (Fernandez et al., 2004). Thus, voltage-dependent 
modulation of the selectivity filter/inactivation gate 
by an aromatic residue or hydrophobic residue in S6 
appears to apply to many, but not all, Kv channels. 

The mechanisms of KCNQ1 K + channel inactivation 
are not well understood, but a functional analysis of 
site-directed mutant channels and homology modeling 
approach was used previously to propose a tripartite 
model of inactivation in these channels. Specific mu- 
tations of Leu273 (S5), Val310 (pore helix), or Phe340 
(S6) in KCNQ1 were found to greatly enhance intrin- 
sic inactivation as well as alter channel deactivation 
(Seebohm et al., 2005). V310Aand V310G KCNQ1 chan- 
nels exhibited the most extensive inactivation compared 
with WT channels, whereas L273F also inactivated sub- 
stantially, but at a much slower rate than the other two 
mutant channels. Further study discovered that L273F- 
induced inactivation differed from inactivation of WT 
channels, and molecular dynamic simulations suggested 
that slow inactivation of L273F channels may be caused 
by selective alteration of the conformation of selectivity 
filter site si (Gibor et al., 2007) . Sequence alignment in- 
dicates that Val310 and Phe340 in KCNQ1 are equiva- 
lent to Leu434 and Tyr464 in hEAGl. Leu273 in 
KCNQ1 is located two helical turns above the position 
equivalent to Phe359 in hEAGl; in both channels, 
these S5 residues make direct contact with the key 
pore helix residue (i.e., Val310 in KCNQ1 or Leu434 
in hEAGl). Thus, in two unrelated Kv channels, 
hEAGl and KCNQ1, selectivity filter-mediated inacti- 
vation can be manipulated by mutation of three resi- 
dues located in equivalent or similar positions in the 
channel subunit. 

Summary 

Human EAG1 channels inactivate in a slow voltage- 
dependent manner. The resulting changes in current 
magnitude are relatively minor and do not have any 
known physiologically relevant consequences. However, 
inactivation can be dramatically enhanced by ICA, a 
drug previously shown to have the opposite activity on 
hERGl channels, or by mutation of Tyr464 located in 
the S6 segment at a position known to be a key determi- 
nant of inactivation in several unrelated K + channels. 

Our experimental data and homology modeling sug- 
gest a minimal tripartite molecular model for hEAGl 



inactivation involving interactions between Tyr464 
in S6, Leu434 in the pore helix, and Phe359 in S5. Res- 
cue of hEAGl from Y464A-enhanced inactivation by 
mutating either Leu434 or Phe359 to smaller residues 
suggests that Tyr464 normally maintains the selectivity 
filter in a conformation that favors ion conductance 
and disfavors inactivation. We propose that interaction 
between Tyr464 and Phe359 places a "brake" on inacti- 
vation by limiting Phe359-Leu434 interaction. Muta- 
tion of Phe359 to Ala or Leu alone does not alter 
inactivation gating, but these mutations can prevent 
Y464A-enhanced inactivation. Phe359 and Tyr464 are 
also key determinants of ICA-mediated enhancement 
of inactivation. When mutated to Leu (F359L), ICA 
activates hEAGl, whereas it inhibits F359L/Y464A 
channel currents. Further studies will be required to 
determine if modulation of inactivation by Y464A and 
ICA occurs by convergent molecular mechanisms, and 
whether the molecular model for selectivity filter-medi- 
ated inactivation of hEAGl and KCNQ1 is applicable to 
other Kv channels. 
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